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We study, theoretically and experimentally, disorder-induced resonances in randomly-layered sam- 
ples, and develop an algorithm for the detection and characterization of the effective cavities that 
give rise to these resonances. This algorithm enables us to find the eigen-frequencies and pinpoint 
the locations of the resonant cavities that appear in individual realizations of random samples, for 
arbitrary distributions of the widths and refractive indices of the layers. Each cavity is formed in 
a region whose size is a few localization lengths. Its eigen-frequency is independent of the location 
inside the sample, and does not change if the total length of the sample is increased by, for example, 
adding more scatterers on the sides. We show that the total number of cavities, JV cav , and reso- 
nances, iV res , per unit frequency interval is uniquely determined by the size of the disordered system 
and is independent of the strength of the disorder. In an active, amplifying medium, part of the cav- 
ities may host lasing modes whose number is less than N ICS . The ensemble of lasing cavities behaves 
as distributed feedback lasers, provided that the gain of the medium exceeds the lasing threshold, 
which is specific for each cavity. We present the results of experiments carried out with single-mode 
optical fibers with gain and randomly- located resonant Bragg reflectors (periodic gratings). When 
the fiber was illuminated by a pumping laser with an intensity high enough to overcome the lasing 
threshold, the resonances revealed themselves by peaks in the emission spectrum. Our experimental 
results are in a good agreement with the theory presented here. 



I. INTRODUCTION 



Anderson localization of waves in one-dimensional ran- 
dom media is associated with an exponential decrease 
of the wave amplitude inside the disordered, locally- 
transparent sample, which results in an exponentially 
small typical transmittance T typ cc exp (— L/l\ oc ) <C 1 
(where L is the length of the sample, and h oc is the 
localization length). Another manifestation of Ander- 
son localization is the existence of resonant frequencies, 
where the transmission increases drastically, sometimes 
up to unity. These frequencies correspond to the quasi- 
localized eigenstates (modes, or disorder-induced reso- 
nances) characterized by a high concentration of energy 
in randomly-located points inside the system. 

Even though one-dimensional (ID) localization has 
been intensively studied during the last few decades^ (see 
also^ and references therein), most of the analytical re- 
sults were obtained for mean quantities, i.e., for values 
averaged over ensembles of random realizations. These 
results are physically meaningful for the self-averaging 
Lyapunov exponent (inverse localization length), which 
becomes non-random in the macroscopic limits. For non- 



self- averaging quantities (field amplitude and phase, in- 
tensity, transmission and refection coefficients, etc.) a 
random system of any size is always mesoscopic, and 
therefore, mean values have little to do with measure- 
ments at individual samples. This is most pronounced 
when it comes to the disorder-induced resonances whose 
parameters are extremely volatile and strongly fluctuate 
from realization to realization^. In particular, ensemble 
averaging wipes out all information about the frequencies 
and locations of individual localized states within a par- 
ticular sample, even though these data are essential for 
applications based on harnessing micro and nano cavities 
with high Q-factors. 



A. Disorder-induced cavities and resonances 

Nowadays, photonic crystals are believed to be the 
most suitable platforms for the creation and integration 
of optical resonators into optical networks. To create an 
effective resonant cavity that supports localized high-Q 
modes in a photonic crystal (PC), it is necessary to break 
periodicity, i.e. to introduce a defect in a regular system. 
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As fluctuations of the dielectric and geometrical parame- 
ters are inevitably present in any manufactured periodic 
sample, it could create a serious obstacle in the efficient 
practical use of PCs. Therefore considerable efforts of 
researches and manufacturers go into the control of fluc- 
tuations. Alternatively, if rather than combating imper- 
fections of periodicity, one fabricated highly-disordered 
samples, they could be equally well harnessed, for ex- 
ample, for creating tunable resonant elements. This is 
because ID random configurations have a unique band 
structure that for some applications have obvious advan- 
tages over those of PC£. Contrary to periodic systems, 
resonant cavities inherently exist in any (long enough) 
disordered sample, or can be easily created by introduc- 
ing a non-random element (for example, homogeneous 
segment) into an otherwise random configuration. The 
effective wave parameters of these cavities are very sensi- 
tive to the fine structure of each sample, and can be easily 
tuned either by slightly varying the refractive index in a 
small area inside the sample or, for example, by changing 
the ratio between the coupling strength and absorption. 
This makes possible to shift the resonant frequency (and 
thereby to lock and unlock the flow of radiation) ; to cou- 
ple modes and create quasi-extended states, etoZ^. 

It has been shown^ that each localized state at a 
frequency u> = w rcs could be associated with an effec- 
tive, high Q-factor resonance cavity comprised of almost 
transparent (for this frequency) segment bounded by 
essentially non-transparent regions (effective barriers). 
Wave tunneling through such a system can be treated 
as a particular case of the general problem of the trans- 
mission through an open resonator—. The distinguish- 
ing feature of a disorder-induced cavity is that it has no 
regular walls (the medium is locally transparent in each 
point), and high reflectivity of the confining barriers is 
due to the Anderson localization. Moreover, different 
segments of the sample are transparent for different fre- 
quencies, i.e., each localized mode is associated with its 
own resonator. 



B. Random lasers 

If the medium inside such a cavity is amplifying, the 
combination of the optical gain and the interference of 
multiply scattered radiation creates coherent field and 
gives rise to multi-frequency lasing with sharply peaked 
lasing spectrum. Random lasers (RL) are the subject of 
increasing scientific interest due to their unusual prop- 
erties and promising potential applications^—. Unlike 
conventional lasers, where any disorder is detrimental, 
in a random version, scattering plays a positive role in- 
creasing the path length and the dwell time of light in 
the active medium. So far most studies, both experimen- 
tal and theoretical, have concentrated on 3D disordered 
systems and chaotic cavities. One of the grave draw- 
backs of 3D random lasers is their inefficient pumping, 
which is hampered by the scattering of the pumping ra- 



diation in the random medium. A one-dimensional RL, 
which is free from this disadvantage, can be realized ei- 
ther as a random stack of amplifying layers^ or as a set 
of Bragg gratings randomly distributed along a doped 
optical fiber—. In the last case, the wavelength of the 
pumping laser is shifted from the Bragg resonance of the 
gratings, and the fiber is excited homogeneously along 
the whole length of the sample. Both these methods 
reduce noticeably the lasing threshold as compared to 
3D random lasing systems. Because the frequencies of 
the modes and the locations of the effective cavities vary 
from sample to sample randomly, in the most cases they 
are described statistically 2 ^— . However, usually we deal 
with a specific random sample, and it is important to 
know how many modes and at which frequencies can be 
excited in a given frequency range; where these modes 
are localized inside this sample, eto^i. 

Another field of research in which this information 
is crucial has arisen recently after it had been realized 
that Anderson resonances could be used to observe cav- 
ity quantum electrodynamics effects by embedding a sin- 
gle quantized emitter (quantum dot) in a disordered PC 
waveguid o 25 i 26 . In those experiments, the efficiency of 
the interaction between radiation and disorder-induced 
cavities depends strongly on the location of the source 
inside the random sample. Indeed, all QED effects are 
well-pronounced when the emitter with a given frequency 
is placed inside the effective cavity, which is resonant at 
this frequency, and could be completely suppressed oth- 
erwise. 

In this paper, we develop an algorithm that enables to 
detect all cavities and to find their locations and eigen- 
frequencies for any individual sample with given geome- 
try and optical parameters. It is shown that in the case 
of uncorrelated disorder, the number of disorder-induced 
resonances per unit frequency interval is independent of 
the strength of the fluctuations and is uniquely deter- 
mined by the size of the random sample. The results 
have been checked experimentally using RLs based on 
a single-mode fiber with randomly distributed resonant 
Bragg gratings developed mi£. 



II. DISORDER-INDUCED CAVITIES AND 
RESONANCES 

It has been shown ii&l&d! that for a quantitative de- 
scription of the wave propagation through a disordered 
sample, it is advantageous to consider it as a random 
chain of effective regular resonators with given Q-factors 
and coupling coefficients. The typical size of each res- 
onator is of the order of the localization length, l\ oc , and 
their centers arc randomly distributed along the sam- 
ple. For manufacturing RLs and for the ability to tailor 
their properties, it is important to know the location of 
the resonant cavity for each eigen-frequencyii. To this 
end, a criterion is necessary, which enables to determine 
whether a given area of a disordered sample is either a 
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resonant cavity or a strongly reflecting (typical) random 
segment. 

To derive such a criterion for randomly layered media, 
we consider a disordered sample consisting of N homo- 
geneous layers with thicknesses dj and dispcrsionless re- 
fractive indices n,j (j = 1,2,..., N) that are statistically 
independent and uniformly distributed in the intervals 
(do — Sd, dg + Sd) and (no — Sn, no + Sn) , respectively. (The 
effects of correlation between the thicknesses of the ad- 
jacent blocks have been considered in^Zr— ). In a system 
with uncorrelated layers, the interface between the j-th 
and (j + l)-th layers is located at a random point Zj and 
can be characterized by complex transmission, tj, and 
reflection, rj, coefficients, which are also randomly dis- 
tributed in the corresponding intervals. The numeration 
is chosen from left to right so that Zj < Zj+\. It is as- 
sumed thereafter that the optical contrast between layers 
is small. Therefore all Fresnel reflection coefficients rj are 
also small 

N«i (i) 

and, consequently, 1— \tj\ 2 <C 1. For the sake of simplicity 
these coefficients will be assumed real. Generalization to 
complex-valued rj and tj is straightforward. 

In what follows, it is convenient to turn to the optical 
lengths, i.e. to replace the thicknesses dj and the coordi- 
nates Zj by 

= n jdj, 

Zj = '^n l z l , (2) 

i<j 

where nj is the refraction index of the j-th layer. Note 
that in these notations the wave number has the same 
value k = k = oj/c in all layers. 

To distinguish between an effective resonator and a 
Bragg reflector, a proper physical quantity representa- 
tive of distinctive properties of these objects has to be 
found. To do this, let us consider a wave that propagates 
rightward from a point z = Zj + [right side of the in- 
terface between the (j — l)-th and j-th layers], where its 
amplitude is A^q (thereafter superscripts (+) and (— ) 
denote the amplitudes of the waves propagating to the 
right or to the left, respectively). After traveling through 
the j-th layer, the wave is partially reflected back from 
the interface between the j-th and (j + l)-th slabs. The 
transmitted part extends through (j + l)-th layer and 
is partially reflected from the next interface, etc. (see 
Fig. EE). 

If the scattering at the interfaces between layers is 
weak, Eq. ([1]), the amplitude ' of the total field re- 
flected from N s layers and returned back to the point 
z = z j + can be calculated in the single-scattering ap- 
proximation, and is equal to: 

A^ = r f\k)A^\ (3) 
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Figure 1: (color online) Schematic diagram of closed trajec- 
tories in a randomly-layered sample. 

where 

r j +) ( fc ) = r i+™ CX P \- 2ik ( S J+ m ~ ' ( 4 ) 

TCI— 1 

and k = uj/c is the wave number. The amplitude Aj = 
rj \k) A^ Q ' is introduced in the same way, with the left 
reflection coefficient rj \k) 

r'f'ik) = - rj- m exp [2ik {zj - %-m)] ■ (5) 

m— 

It is taken into account in Eq. ([5) that the reflection co- 
efficients n and r r for waves that are incident on the 
same interface from the left or from the right have op- 
posite signs: r; = —r r . The field of all waves that made 
a closed path and returned back after consequent reflec- 
tions from N s layers located on the right, and N s layers 
located on the left from the j-th layer, has an amplitude 
A^[ . In the general case, A^' is not equal to the initial 
amplitude A^\ and the difference is 

5A\ +) = A$ 4' = Ajp [1 - A, (k)} , (6) 

where 

A j (k)=r^(k)r^(k). (7) 

The function Aj(fc) is an important characteristic, which 
uniquely determines the resonant properties of any one- 
dimensional wave system. In particular, the eigen- 
numbers can be found as poles of the Green function, 
i.e. as the roots of the equatio n 30 ' 31 

Aj(k)-l^r^(k)r^(k)-l = 0. (8) 

In the case of a closed resonator 

A(fc) = exp (2i7m) , n = 1,2,3 ... , 

i.e. arg A(k) = 2im, and ReA(fc) = 1. For an open 
resonant cavity, the roots k ICS of Eq. ([5]) are complex, 
fc res = kn + ikj. If the Q-factor of a resonator is large 
enough, then kn 3> fcj and 

ReA(fc fl ) = l-k R £/Q >0, (9) 
ImA(fc fl ) ~ 0, (10) 
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where I is the resonator length. 

Note that equation (|T0|) is fulfilled also when 
argA(fc) = irn, which is the Bragg reflection condi- 
tion. In this case, in contrast to a resonant cavity, 
sgn{r' + )(fc)} = — sgn{r' - )(fc)} and the real part of A(k) 
is negative, ReA(fc) < 0. 

These properties of the quantity A(k) are quite general 
and can be used to characterize randomly layered sys- 
tems, in particular, to detect effective resonant cavities 
inside them. Indeed, when for a segment of 2iV s layers 
centered at a point Zj inside a long (N ^> N s ) disor- 
dered sample the imaginary part of A(fc) is zero at some 
k = lu/c, this area is either a resonant (at the frequency 
to) cavity or a localization-induced resonant Bragg reflec- 
tor. What it is indeed is determined by the sign of the 
real part of A: in a resonator Re A > 0, while for a Bragg 
grating Re A < 0. 

The last condition is easy to understand if we notice 
that rj(k) is (— 2fc)-Fourier harmonics of the function 



m— 1 



and A ' (fc) is the (+2fc)-Fourier harmonics of the func- 



\ r i Wl = \ r j\VNioc ~ 1) an d therefore 



iVi 



loc 



(13) 



The characteristic scale of the variation of (k) along 
the sample measured in number of layers is of the order 
of TVioc » 1. Thus, the coefficients rj- (fc) and hence the 
quantities Aj(k) are smooth random functions of j, i.e. 
of the distance along the sample. As the length ./V of a 
sample is large enough, N 3> -/Vi oc , the number of regions 
where Re Aj(fc) > is approximately equal to the num- 
ber of regions where ReAj(fc) < 0. The characteristic 
size of these areas is the localization length. Therefore 
the expected number of cavities, iV cav , resonant at a given 
wave number k is: 



L 

2^1oc 



(14) 



tion 



In order to estimate the number of resonances, N ies , 
in a given interval Afc of the wave numbers (in a given 
frequency interval Aw = cAk), it is necessary to estimate 
the cavity "width", 5k res , in the fc-domain. To do this, 
we note that the variation 5k of the wave number leads 
to a variation 5tpj of ipj = arg Aj(fc): 



A'., 



F] >(*) = £ i~ r 3-m) S[z - (z 3 - m - %)]• (12) 



m=0 



This means that (fc) and A 



j (fc) are Bragg reflection 
coefficients from N s slabs located to the right and to the 
left from the j-th layer, respectively. 

Since our prime interest here are disorder-induced res- 



onators for which 



.(±) 



1, and because in the local- 
ized regime the reflection coefficient from a typical re- 
gion is close to unity when its length is of the order (and 
larger) of the localization length, for further calculations 
we have to choose N s = l\ oc /do = N\ oc . Important to 
stress that each resonator is formed in an area of the size 
of a few localization lengths and is practically unaffected 
by the outer (to this area) parts of the sample. This 
means that the eigen-frequency of an effective cavity is 
independent of it location inside a sample, and does not 
change if the total length of the sample is increased, for 
example, by adding more scatterers at its edges. When 
the thicknesses of the layers are uncorrelated and the 

mean thickness dj = do is large compared to the wave- 
length, kda 3> 2ir, the localization length is determined 
by the mean value of the local reflection coefficients r^. 

To prove this statement, note that the functions r^{k), 



Eqs. (|4|5[) . are sums of -/Vi oc uncorrelated random com- 
plex numbers, and therefore can be described in terms of 
a random walk on the plane (Rcr, Imr), with the single 

step equal to \rj\. Then, the mean absolute value |r*- (k)\ 
is the mean distance from the origin after -/Vi oc steps: 



5 k 



dk 
(+)/ 



Im 



Im- 



1 dry>{k) 



,(+) 



dk 



Im- 



1 dAj(k) 
Aj(k) dk 

1 



dr\-\k) 



,(-) 



dk 



(15) 



The resonant wave number k r is defined by the condi- 
tion (|10p . therefore the second resonance appears in a 
small vicinity of the same layer j when the variation of 
arg Aj(k) approaches 2ir: 5(fj ~ 27T. It is easy to see that 
the largest contribution to 5(fj comes from the layers that 
are the most distant from the j-th layer: 



T d 



logrj ±} (fc) ~ \zj -z j±Nl J ~ 

Mocdo = hoc- (16) 



Thus, 5ipj(k) ~ 2 5kl\ oc 



and the characteristic interval 
5k res between resonant wave numbers localized around 
an arbitrary point Zj can be estimated as 



5kr, 



Moc 



(17) 



In^ this result has been obtained for modes located 
around the center of the sample. 

Equations (fT4|) and (jTTj) allow estimating the number 
of resonances N rcs in the given frequency interval Au; = 
cAfc in the sample of the length N: 



iV re 



Afc 
~5k 

uteres 



NAkd /2ir. 



(18) 



5 



The number of cavities, N cay , and the spacing between 
the resonances, Sk Tes , are inversely proportional^ to ^ oc , 
i.e., it increases when the strength of the scattering in- 
creases. Therefore, it follows from Eq. (|18p that N les does 
not depend on the localization length. In other words, 
the number of resonances (and therefore the number of 
peaks in the transmission spectrum and the number of 
regions with enhanced intensity) in a given frequency in- 
terval arc proportional to the size of the random system 
and are independent of the strength of disorder. 

An important point is that Eq. (fT8"]) gives the total 
number of disorder-induced resonances existing along the 
whole length of a random system in a given range Afc. 
When random lasing is concerned, each sample is an ac- 
tive medium, and a new parameter - the specific gain 
rate g — should be involved. In evaluating the number of 
lasing modes, this parameter has to be compared with 
(7) 

the lasing threshold g c of each cavity, which is different 
for different ones. As it has been shown in^ and men- 
tioned in the Introduction above, each disorder-induced 
resonator occupying an area %-jv s < z < zj + n s is built 
of strongly reflecting (as a result of Anderson localiza- 
tion) effective barriers that confine an almost transparent 
(for the given resonant wave number k) region. Reflec- 
tion coefficients of the left and right (from the center) 
parts of this structure are large, which means that the 
normalized amplitudes f^>\ of the ±2/c-harmonics of the 

of the scatterers in 



f (±) 

J^2k 



1. However, 



distributions F^ ±] [Eqs. HI]), (12 
these parts of j-th cavity are large 
the value jik of the amplitude of the 2fc-harmonics of the 
distribution Fj(z) = F l - + \z) — Fj(z) of the scatterers 
along the whole cavity, can be small, or even equal to 
zero. Indeed, since 



w), including the resonant ones. These harmonics provide 
distributed feedback, as it occurs, for instance, in conven- 
tional (not random) semiconductor distributed feedback 
(DFB) laser (see, e.g.^). As in a conventional DFB laser 
with regular Bragg grating (regular periodic spatial dis- 
tribution of the scatterers) , in a disorder- induced cavity a 
lasing mode is excited when the gain rate g of the medium 
exceeds the value g c ^ of the threshold of the cavity. The 
difference between these two cases is that the whole peri- 
odic grating is characterized by a single value of g c , while 
different parts of a random configuration have their own 
different thresholds gc , whose values are minimal into 
the cavities. In just the same way as the phase shift be- 
tween two identical halves of the regular grating causes 
to decrease the threshold in conventional DFB lasers^, 

the phase difference between complex amplitudes ^ 

and /I2I reduces the lasing threshold of the 2fc-resonant 
random cavity. From the considerations presented in the 
previous paragraph it follows that the information about 
g^' is also contained in the quantities A^: the greater 

Re Aj (k) is for a given cavity, the smaller is its threshold 

(i) (?) 
g c . Thus, only those cavities whose thresholds g c ' are 

less than the medium gain g (i.e. the values ReAj(k) 
exceed the critical value A c ), contribute to the forma- 
tion of the lasing spectrum. Therefore the number of 
lasing modes (number of lines in the lasing spectrum) is 
definitely smaller than the total number of cavities N Tes . 

To conclude this section we note that approach pre- 
sented above could be used in studying various wave sys- 
tems, for example, THz tunable vortex photonic crystals. 
Indeed, many calculations on this problem have already 
been performed^, including the role of disorder, but this 
is beyond the scope of this paper. 



J2k - J_ 2k 



III. NUMERICAL SIMULATIONS 



If. 



21, 



and 



.(+) 



.(+) 



.(-) 



.(-) 



2ReA(fc), 



the value of I/2&I is small when ReA(fc) is close to unity. 
For given values of |r^ + ^| and j^'j, the amplitude \ f2k\ 
is smaller the greater is ReA(fc). In the extreme case, 
the amplitude \f2k\ turns to zero, i.e., the 2/c-harmonics 
is completely suppressed. This happens when the ampli- 
tudes have equal absolute values but opposite signs, 
i.e. are shifted in phase by tt. Thus, the (positive) value 
of Re A(fc) can be treated as a measure of the amplitude 
of the 2fc-harmonics of the cavity and of the phase shift 
between the 2A:-harmonics of its left and right parts. 

In a real random sample, the 2A:-harmonics of the spa- 
tial distributions F- (z) of scatterers arc not equal to 
zero in any effective cavity and for any k (any frequency 



To verify the results of the previous section, we have 
calculated the function Aj(fc) numerically for different 
random samples, and have found the distribution of the 
areas with ReAj(fc) > and RcA_,(fc) < along each 
sample for different wave numbers. These distributions 
have been mapped on the coordinate - wave-number 
plane in Fig. [2] 

Then, the eigenvalues, fc res , have been determined from 
Eq. ([8]), and spatial distributions of the intensity cre- 
ated inside the samples by the incident waves with the 
corresponding resonant frequencies (w res = cfc res ) have 
been found and compared with the map. To facilitate 
this comparison, we take into account that the inten- 
sity pumped by the incident wave into a cavity depends 
on its location inside the sample: it decreases exponen- 
tially as the distance of the cavity from the input end 
increase d 9 ' 34 . As an example, the intensity distributions 
along the same sample illuminated by the same wave, ei- 
ther from the left (blue curve), or from the right, 
(red curve) ends are shown in Fig. [5^,. 
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Normalized coordinate along the sample 

Figure 2: (color online) Mapping of Aj(fc) and the normal- 
ized resonant intensities I s (j, k) on the coordinate-frequency 
plane. The units in the ir-axis are normalized to the local- 
ization length. In the black areas ReAj(fc) > 0, gray color 
corresponds to ReAj(fc) > 0. Approximately half of the full 
length of any line k = ui/c — constant is occupied by black 
areas where the cavities are located. Yellow color marks the 
regions where |argAj(fc) — 2n\ < 0.05. The regions with 
strong concentration of the wave field, where I s (j,k) > 1/2, 
are marked in green. As can be seen, practically all of them 
are located in the cavities (black areas) as predicted here. In- 
set: enlarged view of the selected area with three [(a), (b), 
and (c)] resonances. 



In both curves shown in Fig. |3l the intensity in the 
closest to the input cavity is far above the intensity inside 
the more distant one. Therefore we have introduced the 
normalized quantity 









_max{/,(j)} max{J r (j)}_ 


max | 


hU) | IrU) 


} 


_max{/j(j)} max{J r (j)}_ 



(19) 



which is symmetric with respect to the direction of in- 
cidence and reveals equally well all resonators, indepen- 
dently on their positions Fig. [3Jd. 

In Fig. O regions where ReAj(k) > and ReAj(k) < 
are marked by black and gray, respectively. The nor- 
malized intensity I S (J, k) is shown by green. Yellow color 
marks the regions where |argAj(fc) — 2ir\ < 0.05. One 
can see that any horizontal cross-section k = constant 
contains approximately the same number of black and 
blue regions, and practically all resonances are located 
in "black" areas associated with cavities, as it was pre- 
dicted. 

Examples of the spatial distributions of Re Aj (blue 
curves) for three resonant frequencies marked by (a), 
(b), and (c) in Fig. [2] are presented in Fig. along 
with the corresponding distributions of the intensity (red 
curves). Shown in Fig. [5] cross-correlation functions, 
C(l) = Re Ajlj-i, of those two types of curves re- 




el 500 1000 1500 2000 

Layer number (Dimensionless coordinate along the sample) 



Figure 3: (color online) (a) Spatial distributions of the inten- 
sity generated by a wave incident from the left (thin blue) and 
right (thick red) sides of the sample; (b) intensity normalized 
in accordance with Eq. (|19[) . The incident wave generates a 
resonance (peak of the intensity) only in the cavity closest to 
the edge it is coming from. The distribution of the normal- 
ized intensity is independent of the direction of incidence and 
clearly reveals all cavities resonant at a given frequency. 



veal strong correlation between the normalized intensity 
and Re A. One can see that the cavities are well detected 
by the Re A > criterion. 

Figure [6] demonstrates that, in accordance with 
Eq. (|18|) . the number of resonances in a given frequency 
interval is independent of the strength of disorder (of the 
values of the local reflection coefficients rj) and of the 
variance of the fluctuations of the thicknesses dj , and is 
completely determined by the size of the sample. 

This rather counter-intuitive result is also supported 
by the numerical simulations presented in Fig. where 
the locations of the resonant cavities in the coordinate- 
frequency plane are shown (marked in black) for two sam- 
ples. The samples are geometrically identical, i.e. have 
the same spatial distribution of the scatterers, and dif- 
fer only in the amplitudes of their reflection coefficients, 
so that the localization length at the upper picture is 
twice larger than in the lower panel. When comparing 
both images in Fig. it is easy to see that in passing 
from one picture to another the sizes of the black areas 
in the x-direction increase while the distances between 
them along the y-axis decrease in the same proportion, 
so that, in keeping with Eq. (|18|) . N Tes remains the same, 
meaning that there is no dependence of the number of 
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Figure 5: (color online) Cross-correlation functions C(l) = 
X]., Re Ajlj-i of the normalized intensity and Re A for three 
resonant wave numbers marked by (a), (b), and (c) in the 
inset in Fig. [2] The curves reveal strong correlation between 
the normalized intensity and Re A. 



Figure 4: (color online) Spatial distributions of ReAj (thin 
blue) and of the normalized intensity (thick red) for three 
resonant wave numbers marked by (a), (b), and (c) in the 
inset in Fig. [2] One can see that each resonance is localized 
in the area where Re Aj is positive and takes maximal value. 



resonances on the localization length. 



IV. EXPERIMENTS WITH ACTIVE RANDOM 
SAMPLES: NUMBER OF LASING MODES 

The experimental detection of disorder-induced cavi- 
ties and resonances is a challenging task in optics. Al- 
though the intensity distribution cannot be measured di- 
rectly, the resonances, in principle, can be revealed by 
transmission experiments at different frequencies. How- 
ever, although resonant transmission in lossless systems 
is essentially higher than at typical frequencies, it is much 
stronger affected by absorption, which is proportional to 
the exponentially large intensity inside resonators. In mi- 
crowave cxpcrimcnts&i, where the absorption length was 



much larger than the total length of the system (single- 
mode waveguide) , the transmission was bellow the noise 
level even at centrally located resonances. This fact was 
of little concern in that case because the microwave probe 
could be inserted in different points inside the waveguide. 
In fiber optic systems, however, this is not possible. With 
long optical fibers, to make the transmission measurable 
it might be necessary to compensate for the absorption 
by, for example, introducing amplification. When the 
amplification is sufficiently large, all resonances manifest 
themselves as sharp lines in the transmission spectrum. 

Experimental studies of resonances in disordered op- 
tical fibers with active elements is also of interest for a 
better understanding of the physics of RL, because in am- 
plifying media all resonant modes are potential candidate 
for lasing. The ability to monitor and to tailor their pa- 
rameters, especially the total number, eigen frequencies, 
and locations is critical to fulfill this potential. 

Random one-dimensional cavities can be created in op- 
tical fibers by the introduction of randomly-positioned 
reflectors, in the form of Bragg gratings^. To create 
such structures, we have used commercial Er/Ge co- 
doped fibers (INO Quebec, QC, Canada) that are single- 
mode at 1535 nm. Erbium is the active element and 
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Figure 6: (color online) Number of resonances iVres as a func- 
tion of the averaged local reflection coefficient (strength of 
disorder) (\r\). It is easy to see that, as predicted by the the- 
ory, the number of resonances in a given frequency interval 
grows linearly with the size of the sample and is independent 
of the strength of the disorder. Parameters of the numerical 
simulations: kdo — 10.0, Afcdo = 0.1. Each point is obtained 
by averaging over 10 3 random samples. 
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Figure 7: Distribution of the resonant cavities in the 
coordinate-frequency plane for two geometrically identical 
samples with different strengths of the scattering and thus 
with different localization lengths. Top: (\rj\) = 0.2; bottom: 
(|Tj|) = 0.15. The distances between the black tick marks un- 
der each picture are equal to the corresponding localization 
lengths. 



germanium doping provides a kind of photosensitivity 
that can be used to change, locally, the refractive index. 
The Bragg gratings were fabricated by exposing the fiber 
to UV light (244 nm) from a frequency doubled argon 
ion laser through a periodic mask whose spatial period 
is 1059. 8nm. Each one of the fabricated gratings had a 
length of approximately 5 mm. The random distances 



between the gratings were statistically independent and 
uniformly distributed in the interval do ± 0.8 mm, where 
the mean distance between gratings do was approxi- 
mately 5 mm. 

The fabricated Bragg grating have a narrow reflection 
spectrum (its full width of half maximum is about 0.17 
nm) centered at 1535.3 nm, with a maximum reflectivity 
of about 0.07-0.08. Under these conditions, the estimated 
localization length is found to be about 5-6 gratings. We 
notice that variations in the mask alignment, recording 
exposure, and fiber tension during the writing process 
caused small variations in the central wavelength of the 
gratings and on the sharpness of their reflection spectra. 
As a result, the half width of the reflection spectrum of 
an array of 31 gratings is about 0.27 nm. 

The optical arrangement employed to fabricate the 
gratings and to measure the reflection spectra of the ar- 
rays is illustrated in Fig. [5] Laser action was obtained by 
end-pumping the system with 980 nm radiation from a 
semiconductor laser. A wavelength-division multiplexing 
(WDM) was used to separate the pumping wavelength 
from the radiation emitted by the laser (see Fig [5]) . Mea- 
surements of the gratings transmission/reflection coeffi- 
cients with a spectral resolution of 0.001 nm were carried 
out in the spectral range 1520-1580 nm, using a tunable 
semiconductor laser (New Focus Vclocity-6300) with a 
coherence length of a few meters. As illustrated in the 
figure, new gratings were fabricated in the sequence, be- 
ginning from the pumping end of the fiber. 
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Figure 8: (color online) Schematic diagram of the experimen- 
tal setup. In (a) we illustrate the configuration of the laser. 
The Wavelength-Division Multiplexing (WDM) device per- 
mits the separation of the pump (980 nm) and lasing radi- 
ation (1550 nm). In (b) we show the arrangement used to 
measure the reflection spectrum of the grating array during 
the fabrication process. 

To explore the dependence of the total number of res- 
onances, ./Vrcs, on the size of the system, we measured 
the frequency spectrum of the reflection coefficient at 
samples with different numbers of reflectors, and made 
use of the fact that each resonance manifested itself as 
a sharp drop of the reflectivity. Black points in Fig. [9] 
represent the total number of the resonances detected in 
the arrays of different numbers of gratings N g (different 
lengths of the samples, L = N g do, d$ = 5 mm) in the 
wavelength range (1534.8-1535.6 nm). The theoretical 
prediction Eq. (TT5)) , 7V rcs = 1.67iV s (dotted line), is in an 
excellent agreement with the experimental data. 

In Fig. 1101 we present the results of the measurements 
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Figure 10: (color online) Emission spectra of the samples 
with 7, 14, 20, and 31 gratings for 20 mW (green) and 40 
mW pumping powers. Both values are well above the lasing 
threshold. The vertical lines mark the spectral positions of 
the emission lines. 



Figure 9: Top: Experimental reflection spectrum from an 
array of fourteen gratings. Bottom: Experimental (black 
points) and theoretical (dotted line) dependence of the num- 
ber of resonances, N ICS , on the number of gratings N s . The 
theoretical dependence is calculated using Eq. [18] and the pa- 
rameters of the experiments. 



of the emission spectra of the RL fiber containing 7, 14, 
20, and 31 randomly distributed Bragg gratings (RDBG), 
for two values of the pumping power: 20 mW denoted by 
the continuous line curves, and 40 mW, denoted by the 
dashed line curves. With these pump levels, the systems 
were above threshold in all cases. For the measurements 
we used a spectrum analyzer with the resolution of 0.01 
nm. 

One can see from Fig. [TU] that in the considered fre- 
quency range, the emission spectrum contains several 
peaks that reveal the presence of resonant modes. As 
the number of gratings grows so does the number of res- 
onances: there are two peaks for the RL with seven grat- 
ings and seven peaks for the sample with 20 gratings. 
For the RL with seven gratings, the two peaks main- 
tain their positions and relative intensities as the pump 
power increases. For systems with a higher number of 



gratings, the competition between modes produces tem- 
poral fluctuations in the relative strengths of the emission 
lines; these fluctuations also depend on the pump power. 
While the relative intensity of the peaks can change with 
the pump power, their spectral positions (indicated by 
the vertical lines in the figure) remain fixed. 

Another interesting feature that can be observed in 
Fig. [TU]is that, once an emission line appears in a system 
with a low number of gratings, it is likely to reappear 
in a system with a higher number of gratings. One can 
see, for example, that the emission lines observed with 
the system with seven gratings arc also present in the 
systems with more gratings. 

The curves shown in Fig [TT] represent spectra of the 
radiation emitted by a random laser with 22 gratings, 
measured in 1 second intervals. Even at constant pump- 
ing power, the number of well-defined lasing modes and 
their emission intensity fluctuate in time. The emission 
frequencies, however, remain fixed; one can see that they 
always coincide with one of the vertical lines of the figure. 
The observed fluctuations of the intensity of the spec- 
tral lines could not be caused by relatively small (~5%) 
fluctuations of the intensity of the pumping laser, and 
apparently were associated with non-linear effects. In- 
deed, despite the relatively low power of the emission, 
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the field inside the high-Q cavities can be strong enough 
(due to resonance) to generate a Kerr-type nonlinear- 
ity. These effects arc of importance, for example, in dis- 
tributed Bragg reflector fiber lasers that are longer than 
20cm 37 i 38 
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Figure 11: (color online) Emission spectra of the sample with 
22 gratings generated by the fixed pumping power (40 mW) 
at different moments of time. 



The black points in Fig. [12] denote the number of las- 
ing modes measured in the wavelength range A ± AA = 
1535. 3±0. 3 nm in samples with different numbers of grat- 
ings N g (and different lengths of the system; L = N g do, 
do = 5 mm). These numbers were obtained by adding all 
the emission lines present in the emission spectra over an 
extended period of time. The amount of the lasing modes 
growths linearly with the length of the system, although, 
in accordance with the theoretical reasoning above, it is 
always less than the total number of the effective reso- 
nant cavities in the sample (Fig. [9]). 



V. SUMMARY 

To conclude, eigen-modes (resonances) of a randomly 
layered long (L 3> hoc) sample are localized in disorder- 
induced effective cavities of the size of the order of the 
localization length that are randomly distributed along 
the sample. An algorithm for finding this distribution 
in an individual configuration with arbitrary (random) 
parameters was developed based on the calculation of the 
function Aj(k) in Eq. ([7]). It was shown that a cavity of 
an effective size of a few localization lengths, in which a 
mode with k = k Tes can be localized, exists around j-th 
layer when Im Aj(fc rcs ) = and Re Aj(k Te s) > 0. In the 
case of uncorrelated disorder and weak scattering, the 
spacing between eigen-levels and the number of cavities 
iV r es in a given frequency interval does not depend on the 
strength of disorder, and are uniquely determined by the 
size of the sample. The frequency of each resonance is 
independent of the coordinate of the effective cavity, in 
which it is located. The number of lasing modes depends 
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Figure 12: Number of resonances measured in the wavelength 
range A ± A A = 1.5 • 1(T 4 ± 0.85 • 10~ 8 cm (black points) 
as a function of the number of Bragg reflectors. Since the 
total length of the sample is proportional to the number of 
gratings, the figure actually presents the dependence on size 
of the random system. 



also on the ratio between the threshold values g c of 
the individual cavities and the gain g of the medium, 
and is less than N IOB . The theoretical predictions and 
numerical results are in reasonable agreement with the 
experimental data obtained by measuring the emission 
spectra of the random laser based on the single-mode 
fiber with randomly distributed Bragg gratings. 
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